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ABSTRACT: 

One  and  one-half  hour  isentropic  trajectories  obtained  from  an  objective 
computer  technique  for  the  308K,  312K,  320K  and  332K  isentropic  surfaces  are 
investigated  for  the  May  28,  1967,  severe  local  storm  case  in  the  National 
Severe  Storm  Laboratory  (NSSL)  network.    In  conjunction  with  the  trajectories, 
air  parcel  properties  of  vertical  motion,  acceleration,  divergence,  vorticity, 
static  stability  and  total  energy  have  been  studied  with  respect  to  development, 
growth,  movement  and  decay  of  the  storm  system.    Basic  assumptions  are  that 
the  flow  is  frictionless  and  adiabatic  over  the  entire  grid.    The  results  show 
that  the  observed  winds  do  provide  some  skill  in  formulating  isentropic  trajec- 
tories which  demonstrate  the  presence  of  storm  activity  in  the  specific  region. 
The  investigation  of  a  mesoscale  inertial -gravity  wave  having  a  wave  length 
of  150  KM,  period  of  3  hours  and  phase  speed  of  14  m  sec"1  is  also  presented, 
with  its  possible  correlation  to  the  onset  of  the  severe  storm. 
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ABSTRACT 

One  and  one-half  hour  isentropic  trajectories  obtained  from  an 
objective  computer  technique  for  the  308K,  312K,  320K  and  332K  isen- 
tropic surfaces  are  investigated  for  the  May  28,  1967  severe  local 
storm  case  in  the  National  Severe  Storms  Laboratory  (NSSL)  network. 
In  conjunction  with  the  trajectories,  air  parcel  properties  of 
vertical  motion,  acceleration,  divergence,  vorticity,  static  sta- 
bility and  total  energy  have  been  studied  with  respect  to  develop- 
ment, growth,  movement  and  decay  of  the  storm  system.   Basic 
assumptions  are  that  the  flow  is  frictionless  and  adiabatic  over 
the  entire  grid.   The  results  show  that  the  observed  winds  do 
provide  some  skill  in  formulating  isentropic  trajectories  which 
demonstrate  the  presence  of  storm  activity  in  the  specific  region. 
The  investigation  of  a  mesoscale  inertia 1-gravity  wave  having  a 
wave  length  of  150  KM,  period  of  3  hours  and  phase  speed  of  14  m 
sec  "  is  also  presented,  with  its  possible  correlation  to  the 
onset  of  the  severe  storm. 
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1.   INTRODUCTION 

The  severe  local  storm  is  an  extremely  complex  meteorological 
phenomenon  which  can  lend  itself  to  interesting  and  challenging 
research.   The  approaches  one  can  take  for  extensive  studies  are 
almost  unlimited  as  evidenced  by  the  wealth  of  literature  on  various 
aspects  of  the  severe  storm.   Yet  there  are  many  questions  on  the 
formation,  movement,  growth  and  decay  of  these  systems  which  remain 
unanswered.   Similar  synoptic  conditions  may  seem  to  exist  during 
two  separate  time  periods,  and  yet  storms  may  be  forecasted  to 
develop  and  then  not  form,  or  vice-versa.   In  particular,  the  extent 
to  which  an  individual  storm  may  evolve  is  not  reasonably  understood 
or  predictable. 

A.   SEVERE  STORM  RESEARCH 

In  19^7  a  comprehensive  study  was  undertaken  in  the  Thunder- 
storm Project  which  was  perhaps  the  beginning  of  extensive  scientific 
efforts  into  the  study  of  severe  local  storms.   From  this  project 
Byers  and  Braham  (19^-9)  provided  new  concepts  on  the  structural  fea- 
tures of  such  storms.   The  National  Severe  Storms  Project,  initiated 
in  196l,  was  an  effort  designed  to  study  Great  Plains  convective 
storms  and  the  environment  in  which  they  form.   More  recently  the 
Barbadoes  Experiment,  which  has  provided  data  for  study  of  tropical 
cumulus  convection,  and  Project  Storm  Fury,  for  the  study  of  hurri- 
canes, have  enabled  investigators  to  carry  on  detailed  research  into 
severe  storm  characteristics. 

Mesoscale  networks  located  in  various  regions,  with  instrument- 
ation positioned  at  specified  spacing,  have  provided  for  collection 


of  data  for  certain  research  objectives.   In  particular,  the  National 
Severe  Storms  Laboratory  (NSSL),  Norman,  Oklahoma,  has  a  systematic 
network  of  surface  and  upper  air  reporting  stations  and  radar  facili- 
ties. 

There  have  been  several  studies  involving  kinematic  properties 
which  have  given  insights  into  storm  characteristics.   Williams  (1962) 
investigated  small-scale  kinematic  properties  of  surface  horizontal 
velocity  divergence,  horizontal  stretching  deformation,  and  the 
vertical  component  of  horizontal  shearing  deformation  for  isolated 
thunderstorms.   He  found  the  intensity  and  rate  of  change  of  these 
properties  for  small  scale  systems  were  2  to  3  orders  of  magnitude 
greater  than  those  normally  measured  on  the  gross  synoptic  scale. 
Endlich  and  Mancuso  (1968)  investigated  a  number  of  kinematic 
quantities  using  objective  analysis  procedures  and  compared  them 
with  storm  developments.   They  found  vorticity  production  in  low 
levels  and  horizontal  divergence  in  the  upper  troposphere  to  be 
good  objective  indicators  of  severe  storm  activity.   Fankhauser 
(196M  used  a  statistical  and  observational  approach  to  propose 
that  the  translation  of  convective  complexes  is  related  in  a 
systematic  way  to  the  mean  wind  through  the  tropospheric  layer  in 
which  the  clouds  are  imbedded,  and  that  this  relationship  is  depen- 
dent on  the  unique  character  of  the  mean  wind  field.   Hammond  (1967) 
investigated  an  isolated  airmass  storm,  which  moved  rapidly  north- 
eastward across  central  Oklahoma  where  its  movement  was  to  the  left 
of  all  winds  in  the  middle  and  upper  troposphere.   He  concluded  that 
the  same  environmental  wind  field  can  account  for  the  characteristic 
organization  and  structure  of  either  right  or  left  moving  storms, 
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depending  on  storm  velocity.   Studies  by  Browning  and  Fujita  (1965), 
Newton  and  Fankhauser  (1964)  and  Browning  and  Ludlum  (1962)  have 
shown  that  large  thunderstorms  normally  travel  to  the  right  of  the 
winds  in  the  lower  and  middle  level  of  the  troposphere.   Matsumoto, 
Ninomiya  and  Akiyama  (1967)  studied  the  relations  of  cumulus  con- 
vections and  the  mesoscale  convergence  field  over  the  coastal  area 

of  Japan  and  found  that  well  developed  cumuli  exist  only  in  the  area 

-4    -1 
of  surface  convergence  of  more  than  10   sec 

B.   ISENTROPIC  ANALYSIS  AND  TRAJECTORY  STUDIES 

The  use  of  isentropic  analysis  to  study  the  physical  processes 
at  work  in  the  atmosphere  is  well  known  since  Rossby  (1937)  intro- 
duced the  technique.   The  analysis  and  its  uses  are  well  described 
by  Namias  (1940).   Tongues  of  dry  and  moist  air  from  specific  humidity 
analysis  on  isentropic  surfaces  show  moisture  sources  and  instability 
regions  for  thunderstorm  activity. 

Calculations  of  isentropic  trajectories  for  the  study  of  moisture 
sources,  vertical  motions,  velocity  accelerations,  type  of  air  and 
circulations  within  systems,  vorticity  sources,  and  for  investigations 
of  transports  of  pollutants  and  trace  substances  has  long  been  appeal- 
ing since  non-precipitating  motions  in  the  free  atmosphere  may  be 
assumed  to  be  isentropic  for  limited  periods  of  time.   Many  investi- 
gators have  used  graphical  calculations  from  isotach  and  streamline 
analysis  to  determine  the  trajectories  (Danielsen,  1959;  Reiter 
and  Danielsen,  1960;  Reiter,  1963).   The  use  of  isentropic 
trajectories  in  thunderstorm  studies  (Bradley  and  Feteris ,  1966; 
Feteris  and  Stout,  1966;  Carlson  and  Ludlum,  1968)  can  provide 
useful  information  in  the  study  of  these  mesoscale  systems. 
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It  is  intended  in  this  research  to  investigate,  by  the  use  of 
isentropic  trajectories,  meteorological  properties  of  air  parcels 
and  their  relationship  to  severe  local  storms.  Several  questions 
will  be  considered: 

1.  Do  isentropic  trajectories  over  a  short  period  of  time 
describe  flow  patterns  which  adequately  show  the  formation,  growth, 
movement ,  and  decay  of  the  severe  storm? 

2.  Do  associated  parameters  computed  along  isentropic 
trajectories  describe  motions  directly  related  to  the  dynamics  of 
the  severe  storm? 

3.  Does  a  combination  of  any  computed  parameters  correlate 
well  with  the  storm? 

h.      Is  there  an  oscillation  or  wave  phenomenon,  in  particular 
a  gravity-inertial  wave,  that  can  be  detected  from  the  data  and/or 
computed  parameters  which  might  relate  to  the  local  severe  storm? 

It  is  intended  to  deal  with  these  specific  questions  in  the 
ensuing  sections  by  conducting  a  case  study  of  a  specific  severe 
storm  system.   The  isentropic  trajectories  will  be  computed  from 
an  objective  computer  technique  in  order  that  a  sufficient  number 
of  trajectories  and  isentropic  surfaces  can  be  studied  for  the  severe 
storm  case. 
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2.   OBJECTIVE  ISENTROPIC  TRAJECTORY  COMPUTER  MODEL 

Danielsen  (1961)  introduced  an  objective  method  for  computing 
isentropic  trajectories  based  on  the  total  energy  equation.   This 
provided  the  foundation  for  the  technique  for  which  Mahlman  (1968) 
developed  a  computer  method  for  an  objective  computation  of  adia- 
batic  three-dimensional  trajectories. 
A.   BASIC  DESIGN 

Utilizing  the  graphical  method  by  Danielsen  (1961),  the 
computer  program  applies  an  objective  method  which  confines  the 
computed  trajectory  from  deviating  too  far  laterally  from  its  true 
position.   The  concept  of  the  system  eminates  from  the  horizontal 
equation  of  motion  in  isentropic  coordinates,  which  is 

^  =  -  f  Ik  x  \V  -  VQM  +IF  (1) 

where  \V  is  the  two  dimensional  horizontal  wind  vector,  f  the 
Coriolis  parameter,  |k  the  unit  vector  in  the  vertical,  V  the  two 
dimensional  gradient  operator  on  a  surface  of  constant  potential 
temperature  0,  M  the  Montgomery  stream  function,  and  |F  the  friction 
force. 

The  basic  model  utilizes  variables  of  true  wind  direction,  wind 
speed,  pressure  and  Montgomery  stream  function  at  each  point  in  a 
four-dimensional  (x,  y,  9,  t)  grid.   The  first  step  in  the  program 
is  to  obtain  u  and  v  components  of  the  velocity  field,  with  the 
effects  of  latitude/longitude  type  griding  included  in  the  calcu- 
lations over  the  data  area.   Then,  the  u  and  v  fields  are  linearly 
interpolated  to  give  values  at  twelve  time  step  intervals. 


13 


The  first  approximation  of  the  trajectory  is  using  a  straight- 
forward kinematic  approach.   Beginning  at  a  grid  point,  a  parcel  is 
displaced  for  one  time  step  using  the  initial  values  of  u  and  v  to 
determine  the  length  and  direction  of  the  displacement.   This  dis- 
placement is  determined  by: 

S  =  \VAt  (2) 

where  S  is  the  displacement  and  \V  is  the  average  wind  velocity 
along  the  trajectory  over  a  At  interval.   Then  the  time  interpolated 
u  and  v  values  are  spatially  interpolated  from  the  four  nearest  grid 
points  to  the  parcel  to  obtain  new  u  and  v  wind  components  at  the 
parcel's  new  position.   These  new  values  are  then  used  to  extrapolate 
the  parcel  for  the  next  time  step.   This  process  is  repeated  until 
twleve  time  steps  are  reached.   At  each  point  along  the  trajectory  a 
test  is  made  to  determine  if  the  trajectory  remains  within  the  original 
grid.   If  not,  then  the  trajectory  is  terminated  at  the  point  where  it 
leaves  the  field  and  a  new  trajectory  is  begun  at  the  next  grid  point. 
These  isentropic  kinematic  trajectories  can  provide  a  valid  result  as 
long  as  the  original  wind  fields  have  been  adequately  specified,  and 
the  time  interpolations  are  satisfactory. 

It  is  important  to  note  that  for  inclusion  of  changes  in  the 
Montgomery  stream  function,  Mahlman  (1968)  indicated  that  the  cor- 
rections were  sensitive  to  analysis  errors  in  the  M  field,  rapid 
changes  of  wind  direction  perpendicular  to  the  estimated  trajectory 
path,  and  for  cases  where  the  wind  speed  is  less  than  10  m  sec-  . 
Since  preliminary  studies  of  the  data  indicated  cases  where  the  wind 
flow  was  perpendicular  to  isolines  of  Montgomery  stream  function  and 
wind  speeds  less  than  10  m  sec"  ,  the  subsequent  research  was  thus 
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conducted  from  the  isentropic  kinematic  trajectories.   The  incorpo- 
ration of  adjustments  for  changes  in  Montgomery  stream  function  is 
a  separate  subroutine  in  the  program  so,  if  only  kinematic  tra- 
jectories are  desired,  the  subroutine  can  be  left  out.   The  program- 
ming is  included  in  Appendix  B. 

B.   MODIFICATIONS  FOR  SEVERE  STORM  RESEARCH 

An  additional  field  of  relative  humidity  was  added  in  order 
to  include  a  moisture  constituent.   Routines  to  compute  certain 
meteorological  parameters  were  incorporated  so  that  initial  con- 
ditions and  variations  along  the  kinematic  isentropic  trajectories 
could  be  studied.   These  routines  include  calculations  of: 


_  fev  _  3u\ 
£x"  "  dy) 

3 


+  f  (3) 


9    W        dy) 

A  =  £u_^  (5) 

Bfi  =  §7+f^  <6> 

e    ^x   dy)  Q 

39  (7) 

o  =    -   — 

3p 
where  r   is  the  absolute  vorticity,  Dfl  the  horizontal  divergence, 
Aq  the  stretching  defor;  ation,  Bfi  the  shearing  deformation,  and  a 
the  static  stability.   Also  the  height  (z),  specific  humidity  (q), 
saturated  specific  humidity  (q.s)  and  temperature  (T)  fields  for 
each  isentropic  surface  and  data  time  were  determined.   Along  each 
kinematic  trajectory  the  initial  and  final  values  of  specific  humidity, 
saturated  specific  humidity,  divergence,  absolute  vorticity,  Montgomery 
stream  function,  kinetic  energy,  total  energy,  static  stability,  and 
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7 
potential  vorticity  were  computed.   Total  energy,  in  units  of  10 

ergs  gm   ,  was  defined  as: 

v2 

TE  =  C  T  +  gZ  +  Lq  +  j—  (8) 

P  2 

where  C   is  the  specific  heat  at  constant  pressure  and  L  the  latent 
P 

heat.   Potential  vorticity  was  computed  as: 

P  =  ace  (9) 

where  ^  is  in  units  of  10    K  sec  cm  gm 

The  latitude  and  longitude  positions,  u  and  v  components,  and 
interpolated  height  at  each  of  the  twelve  time  steps  of  the  tra- 
jectory were  obtained.   The  average  vertical  motion  (id)  along  the 
trajectory  was  calculated  using  the  procedure: 

Pf  "  P- 

«  =  ^l^  (") 

with  oo  in  mb  hr   and  where  the  i  and  f  refer  to  initial  and  final 

pressure  over  the  interval  At.   The  acceleration  of  the  parcel  over 

the  trajectory  was  also  calculated.   Finally,  the  pressure  at  which 

the  air  parcel  would  become  saturated  (P  )  and  the  time  of  saturation 

s 

(t  )  were  determined  by: 


-*iU°*  cVr8^  (cV-i) 


si       p 
P  -  P. 


(11) 


t  .-S_±At  (12) 

s   ff       ri 

where  the  i  and  f  refer  to  initial  and  final  points  on  the  trajectory 
and  P,  q  and  q  are  as  previously  defined. 

A  routine  to  utilize  the  CALCOMP  plotter  with  the  IBM  OS/360  com- 
puter was  included  in  order  to  graphically  display  the  computed  tra- 
jectories on  ax,  y  grid  system.   The  data  generated  by  the  program 
was  displayed  such  that  subsequent  studies  could  be  made. 
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3.   DATA  AND  STORM  CHARACTERISTICS 

Data  that  was  considered  in  the  case  study  are  from  observations 
in  the  NSSL  mesoscale  network  on  May  28,  19^7-   Upper  air  data  from 
nine  radiosonde  stations  systematically  spaced  over  the  network  region 
(Fig.  1)  provided  data  at  1700,  1830,  2000,  2130  and  2300  GMT. 

Utilizing  computer  processing  of  mesoscale  rawinsonde  data 
(Kreitzberg  and  Brockman,  196l),  fields  of  Montgomery  stream  function 
(10  erg  gm   )  ,  wind  direction  (deg  from  north),  wind  speed  (m  sec-!), 
pressure  (mb)  and  relative  humidity  (%)    from  the  isentropic  plane 
section  were  hand  analyzed  over  the  network  area.   Allowances  for 
balloon  drift  and  non-simultaneous  releases  were  not  taken  into 
account.   However,  a  comparison  was  made  between  the  original  analysis 
and  a  re-analysis  taking  into  account  balloon  drift  and  storm  location. 
These  results  are  reported  on  in  later  sections.   A  9  x  11  latitude- 
longitude  grid  (Fig.  l),  with  grid  separation  of  .25  degrees,  was 
established  (southeastern  boundary  -3*+N,  97W).   Analyzed  data  were 
then  interpolated  to  the  nearest  grid  point.   Theta  surfaces  308K, 
312K,  320K  and  332K  were  those  selected  for  the  analysis  levels. 
The  308K  theta  surface  was  the  lowest  level  for  which  data  were 
available  at  most  stations  over  the  time  periods.   The  332K  theta 
surface  was  selected  as  being  representative  of  the  upper  level  of 
storm  influence.   This  surface  was  located  at  a  height  of  about  300 
mb. 

A.   BRIEF  SYNOPTIC  ANALYSIS 

The  surface  synoptic  situation  over  Oklahoma  at  1800  GMT  on 
May  28,  1967  indicated  a  weak  high  pressure  (1018  mb)  over  the 
Oklahoma  Panhaldle ,  with  low  pressure  areas  located  over  NW  Missouri 
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CHK  -  Chickasha,  Oklahoma 
COR  -  Cordell/  Oklahoma 
FSI  -  Ft,  Sill,  Oklahoma 
LT3  -  Altus,  Oklahoma 
PVY  -  Pauls  Valley,  Oklahoma 
RIN  -  Ringling,  Oklahoma 
SPS  -  Sheppard  AFB,  Texas 
TIK  -  Tinker  hFB,  Oklahoma 
WAT  -  Watonga,  Oklahoma 


Figure  1,   NSSL  uppar  air  network  and  latitude-longtitude 
grid  system,   Stippled  area  represents  storm  region, 
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(1010  mb)  and  in  central  New  Mexico  (1005  mb).   A  weak  cold  front  lay- 
between  the  two  low  areas,  passing  through  central  Oklahoma.   The 
front  remained  nearly  stationary  during  the  period  under  consideration, 
but  began  to  dissipate  and  become  a  warm  front  by  2^00  GMT. 

The  1800  and  2U00  GMT  500  mb  height  synoptic  analyses  showed 
straight  southwesterly  flow  of  about  15  kt.   The  300  mb  height  chart 
indicated  westerly  flow  of  20-30  kt  at  1200  GMT  with  winds  shifting 
to  southwesterly  and  increasing  slightly  by  2^00  GMT. 

As  reported  by  Fankhauser  (1968),  the  characteristics  were  only 
marginally  conducive  to  thunderstorm  formation,  with  slightly  negative 
stability  indexes  on  morning  soundings.   The  stationary  front  provided 
sufficient  low-level  convergence  to  set  off  weak  convective  storms  over 
the  northwest  portion  of  the  network  by  midday. 

B.   STORM  CHARACTERISTICS 

The  initial  radar  presentation  available  for  study  was  at  1817 
GMT.   Small  echoes  were  present  in  the  western  grid  area  about  30-50 
nautical  miles  from  station  TIK  at  this  time.   Larger  and  more  intense 
echoes  were  present  by  2000  GMT.   Storm  tops  to  the  northwest  of 
station  TIK  exceeded  11  KM  in  height,  and  the  storm  system  was  oriented 
SW-NE  across  the  grid  area.   The  2100  GMT  radar  reflection  indicated 
the  system  to  the  west  of  station  TIK  had  decreased  in  intensity  and 
size,  with  tops  reaching  to  about  9  KM. 

Storm  motions  used  for  calculating  the  relative  trajectories 
were  obtained  from  viewing  surface  radar  presentations  of  similar 
echoes.   These  velocities  were  determined  as: 
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TIME  DIFFERENTIAL 

COURSE 

(deg 

from) 

SPEED  (m  sec"  ) 

1700-1830 

GMT 



0.0  (no  echo) 

1830-2000 

268 

1.5 

2000-2130 

2jk 

2.1 

2130-2300 

272 

k.l 

Fankhauser  (1968)  indicated  that  echo  velocity  was  2^7  (deg 
from)  at  a  speed  of  8  m  sec  .  Also  this  course  and  speed  were 
remarkably  constant  throughout  the  lifetime  of  the  storm.  The 
mean  wind  direction  through  the  atmospheric  layer  for  the  1700- 
2300  GMT  period  was  2^7  (deg  from)  and  had  a  mean  wind  speed  of 
8  m  sec 
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k.      DISCUSSION  OF  RESULTS 

A.   SEVERE  STORM  RELATIONSHIPS  TO  BASIC  DATA  ANALYSES 

There  were  five  basic  fields  which  were  analyzed  on  the  isen- 
tropic  surfaces  in  the  mesoscale  network.   No  objective  smoothing  was 
performed  on  these  fields.   Part  of  the  research  was  to  investigate 
these  basic  fields  and  to  determine  if  they  provide  an  indication  of 
storm  development,  growth,  movement  and  decay. 

The  mesoscale  analysis  over  the  NSSL  network  disclosed  that  a 
region  of  lower  Montgomery  stream  function  values  was  maintained 
near  the  middle  of  the  grid  area  and  to  the  southwest  of  the  storm 
region.   This  prevailing  pattern  was  consistent  on  each  theta  surface. 
Values  were  generally  larger  at  the  initial  data  time  of  1700  GMT  and 
then  decreased  in  value  over  the  entire  grid  as  time  progressed.   How- 
ever, the  predominant  contour  pattern  remained  the  same  (values  on  the 
308K  theta  surface  at  1700  GMT  ranged  from  310. Uo  joules  gm~  to 
310.62;  at  2130  GMT  values  ranged  from  310.17  to  310. 1*5).   The  hori- 
zontal gradient  (V.M)  was  stronger  to  the  SE  of  the  storm  area.   A 

a 

stronger  gradient  existed  at  2300  GMT  on  the  theta  surfaces  than  for 
previous  times.   It  was  noted  that  the  thunderstorm  had  almost  dissi- 
pated by  this  time.   A  dome  of  colder  temperatures  was  associated 
with  the  low  stream  function  values.   This  low  M  value  and  strong 
gradient  region  could  be  associated  with  the  stationary  front  which 
existed  over  central  Oklahoma  and  would  indicate  the  thunderstorm  was 
occurring  northwest  of  the  front.   Figure  2a  is  an  example  of  the 
Montgomery  stream  function  pattern. 
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Figure  2a,  320K  Theta  surface*  2000  GMT  Montgomery  stream 
function  analysis  (add  320,00  to  each  value).  Heavy  solid 
line   represents  region  of  storm  activity,    Values   in  joules 
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conducted  from  the  isentropic  kinematic  trajectories.   The  incorpo- 
ration of  adjustments  for  changes  in  Montgomery  stream  function  is 
a  separate  subroutine  in  the  program  so,  if  only  kinematic  tra- 
jectories are  desired,  the  subroutine  can  be  left  out.   The  program- 
ming is  included  in  Appendix  B. 

B.   MODIFICATIONS  FOR  SEVERE  STORM  RESEARCH 

An  additional  field  of  relative  humidity  was  added  in  order 
to  include  a  moisture  constituent.   Routines  to  compute  certain 
meteorological  parameters  were  incorporated  so  that  initial  con- 
ditions and  variations  along  the  kinematic  isentropic  trajectories 
could  be  studied.   These  routines  include  calculations  of: 


I 


D     =    (9u        3v\  (U) 


A     =     &L-M  (5) 

Ae        ^      ^y]Q  [5) 

Bfl  =  lr  +  r\  (6) 

9  ^x        dyj  Q 

36  \  (7) 

c  -    "  — 
3p 

where  z,      is  the  absolute  vorticity,  DA  the  horizontal  divergence, 

Aq  the  stretching  defor;  ation,  Bfi  the  shearing  deformation,  and  a 

the  static  stability.   Also  the  height  (z),  specific  humidity  (q), 

saturated  specific  humidity  (q.s)  and  temperature  (T)  fields  for 

each  isentropic  surface  and  data  time  were  determined.   Along  each 

kinematic  trajectory  the  initial  and  final  values  of  specific  humidity, 

saturated  specific  humidity,  divergence,  absolute  vorticity,  Montgomery 

/  \ 

stream  function,  kinetic  energy,  total  energy,  static  stability,  and 


Pressure  analyses  on  each  theta  surface  exhibited  sporadic 
patterns.   There  was  no  real  evidence  of  vertical  consistency  with 
regions  of  high  or  low  pressure,  and  no  time  dependency  as  to  move- 
ment or  preservation  of  any  defined  center.   There  were  large  pres- 
sure ranges  in  the  grid  area  on  the  308K  and  312K  theta  surfaces  at 
2130  and  2300  GMT  (7I+8-865  mb  for  308K  and  6U9-T53  mb  for  312K). 
These  pressure  ranges  were  a  result  of  the  decrease  in  height  of  the 
308K  and  312K  surfaces  in  the  SW  corner  of  the  grid  at  2130  GMT  and 
extending  into  the  center  of  the  grid  by  2300  GMT.   Pressure  analyses 
on  the  320K  theta  surface  shows  that  a  strong  gradient  existed  to  the 
SE  of  the  storm.   An  example  of  the  pressure  pattern  is  given  in 
Figure  2b. 

An  analysis  of  relative  humidity  revealed  a  region  in  the  S-SE 
area  of  the  network  which  was  a  minimum  value  area  on  all  theta  sur- 
faces at  all  times  (Fig.  2c).   Since  the  storm  activity  occurred 
between  upper  air  network  stations,  there  was  no  clear  indication  of 
a  region  of  high  relative  humidity.   A  moist  or  dry  tongue  was  not 
evident  from  the  relative  humidity  analyses. 

There  was  a  region  of  higher  wind  speed  in  the  NW  corner  of  the 
grid  area  and  lower  wind  speeds  in  the  SE  corner  (Fig.  2d)  on  all 
theta  surfaces.   This  wind  maximum  propagated  toward  the  east  with 
time  and  increased  in  intensity  until  2300  GMT  when  the  speed  began 
to  decrease. 

Wind  direction  was  generally  anticyclonic  in  the  SE  area  and 
cyclonic  in  the  NW  region  with  a  shear  line  along  the  thunderstorm 
activity  (Fig.  2e).   Time  sections  at  stations  LTS,  FSI,  CHK  and  TIK, 
along  with  the  theta  surface  analyses  of  wind  direction,  indicated 
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Figure  2b,   320k  Theta  surface,  2000  GMT  pressure  (rob)  analysis^ 
Heary  solid  line  represents  region  of  storm  activity, 
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Figure  2c,  320K  Theta  surface,  2000  GM?  relative  humidity 
(%)  analysis,  Heavy  solid  line  represents  region  of  storm 
activity, 
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Figure   2d,      320K  Theta  surface,    2000  GMT  wind  speed   (m  sec''  ) 
analysis,      Heavy  solid  line   represents   region  of   storm  activity, 
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Figum   2a,      320K  Thata  «urfaca,    2000  GMT  wind  direction    (deg« 
from  north)   analy»i«,      Heavy  solid  line   represent*   region  of 
■torn  activity, 
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the  wind  veered  with  height  and  with  time  on  the  308K,  312K  and  320K 
theta  surfaces.   The  time  sections  indicated  the  wind  began  to  back 
with  height  above  the  320K  surface.   This  backing  was  evident  by  1830 
GMT  at  station  LTS,  by  2000  GMT  at  FSI ,  by  2130  GMT  at  CHK  and  by 
2130  GMT  at  TIK.   Wind  speed  at  these  stations  generally  increased 
with  time  and  with  height. 

Data  from  stations  LTS  and  FSI,  which  are  toward  the  SW  of  the 
storm  area,  indicated  the  low  level  wind  speed  increased  until  about 
2015  GMT  and  then  began  to  decrease.   Stations  CHK  and  TIK,  located 
in  the  SE  edge  of  the  storm  region,  show  the  wind  speed  continues  to 
increase  over  the  whole  time  period  (1700-2300  GMT).   Both  the  veering 
with  height  and  speed  increase  have  been  correlated  with  severe  storm 
activity  (Miller,  1967). 

B.   TRAJECTORY  RELATIONSHIPS  TO  THE  SEVERE  STORM 

The  kinematic  isentropic  trajectories  were  computed  for  one  and 
one-half  hour  time  periods  at  each  grid  point  and  on  each  theta  surface. 
For  reference  to  the  trajectories  that  are  discussed,  see  Appendix  A. 

The  1700-1830  GMT  trajectories  clearly  show  veering  with  height 
and  time  in  the  wind  field.   The  308K  and  312K  trajectories  indicate 
anticyclonic  flow  over  the  whole  grid,  whereas  the  320K  and  332K 
trajectories  have  cyclonic  flow  in  the  W-NW  corner  and  anticyclonic 
flow  in  the  S-SE  area.   This  is  the  flow  pattern  just  prior  to  the 
storm  formation,  since  the  thunderstorm  did  not  appear  on  radar  until 
about  1817  GMT. 

The  1830-2000  GMT  trajectories  imply  changes  in  the  wind  field 
that  are  significant.   On  the  308K  surface,  the  trajectories  are 
beginning  to  have  cyclonic  flow  in  the  W-NW,  with  a  shear  line  almost 
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along  the  line  of  thunderstorm  activity.   The  312K  trajectories  specify 
anticyclonic  flow  over  the  whole  area.   The  320K  surface  flow  pattern 
exhibits  cyclonic  flow  along  the  south  edge  of  the  storm  area,  anti- 
cyclonic  flow  to  the  north,  and  straight  flow  along  the  center  line 
of  the  storm  activity.   The  332K  trajectory  pattern  denotes  cyclonic 
flow  to  the  SW  and  along  the  southern  edge  of  storm  activity,  with 
anticyclonic  flow  in  the  NW  and  SE  grid  area.   The  storm  appears  to  be 
oriented  along  the  320K  trajectory  flow  pattern  over  this  1830-2000  GMT 
time  frame. 

During  the  time  interval  2000-2130  GMT,  the  storm  was  in  its  most 
intense  stages.   It  appears  to  be  more  oriented  along  the  308K  flow 
pattern,  which  is  from  the  south  and  cyclonic  over  most  of  the  grid. 
The  area  to  the  west  of  the  storm  is  strongly  cyclonic.   The  312K 
trajectories  indicate  straight  SW  flow.   The  320K  trajectories  give 
evidence  of  a  cyclonic  area  just  to  the  south  of  the  more  intense 
regions  of  the  storm.   The  332K  trajectories  display  anticyclonic  flow 
over  the  entire  grid.   There  is  a  veering  of  70  degrees  or  more  between 
the  308K  and  320K  trajectories  over  this  period  and  pronounced  backing 
in  the  NW  region  between  the  320K  and  332K  surfaces. 

Radar  presentations  indicate  the  dissipation  of  the  storm  over  the 
2130-2300  GMT  period.   The  312K  trajectories  are  from  the  SW  and  show 
slight  anticyclonic  turning.   The  320K  and  332K  trajectories  indicate 
a  slight  cyclonic  turning  in  the  western  portion  of  the  grid  and  anti- 
cyclonic turning  in  the  eastern  portion.   Strong  backing  of  k5   degrees 
or  more  with  height  is  evident  between  the  320K  and  312K  trajectory 
patterns. 

Selected  composite  trajectories  were  obtained  from  each  time  interval 
in  order  to  investigate  air  parcel  flow  over  the  whole  time  period  . 
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These  trajectories  and  associated  parameters  are  discussed  in  the 
following  section. 

C.   RELATIONSHIP  OF  ASSOCIATED  PARAMETERS  TO  THE  SEVERE  STORM 
One  of  the  most  important  parameters  to  be  computed  from  the 
kinematic  isentropic  trajectory  program  was  the  average  vertical  motion 
distribution  for  each  theta  surface  and  each  1.5  hour  time  interval. 
When  used  alone,  these  average  vertical  motion  patterns  are  not  suffi- 
cient to  describe  the  storm  activity.   However,  they  do  provide  very 
definite  clues  to  the  storm  dynamics  when  used  in  conjunction  with 
other  properties. 

During  the  time  interval  1700-1830  GMT,  the  vertical  distribution  of 
vertical  motions  over  the  NSSL  network  indicated  lower  level  air  parcel 
sinking  in  the  western  half  of  the  grid  of  about  8  mb  hr  ""  and  rising 
motions  of  10-20  mb  hr   in  the  eastern  half  of  the  grid.   Upper  level 
patterns  were  reversed,  with  rising  motions  in  the  western  half  and 
sinking  motions  in  the  eastern  half.   In  addition,  there  was  very  strong 
rising  motion  (20-30  mb  hr   )  in  the  western  half  of  the  grid  at  a  height 
of  about  10  KM.   It  is  possible  that  the  high  level  vertical  motion  pat- 
tern is  the  initial  indication  of  an  upper  level  formation  of  supporting 
vertical  motions  for  the  storm. 

The  1830-2000  GMT  vertical  motion  patterns  (Fig.  3a-3d)  clearly 
delineate  and  support  the  storm  system.   At  each  level  there  is  strong 
rising  motion  as  air  parcels  come  into  the  storm  area,  with  sinking 
motions  over  the  remaining  grid  area.   The  320K  isentropic  surface 
vertical  motion  pattern  indicates  the  region  of  maximum  rising  motion. 
The  storm  is  also  oriented  parallel  to  the  320K  trajectory  flow  pattern 
over  this  period.   It  would  appear  the  320K  (i.e.  approximately  500  mb) 
flow  is  responsible  for  maintaining  the  storm  at  this  time. 
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Figure  3a,  308k  Theta  surface,  1830-2000  GMT  ar*rag«  vertical 
motions  (rob  hr"')  and  kinematic  isentropic  trajectories,  Storm 
position  as  of  1915  Gift, 
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Figur*  3b,  312k  That*  mirfaca,  1830*2000  GMT  ar-raga  rertical 
motion*  (lib  hr"4)  and  kimeroatic  i»entropic  tr»j«ctori««,  Stonn 
poaitioa  aa  of  1915  GMT, 
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rigmr*  3c ,  3201  Th»t*  »urf*c«,  1830-2000  GMT  «v»r*g»  vertical 
aotlo««  <ab  hr"*)  &ad  klftaaatlc  i»«*tropic  tr«j«ctori««,  Storm 
poaltloa  *«  of  1915  GMT, 
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Figur*  3d,  332k  Thata  surfaca,  1830-2000  GMT  *Y«ra<ja  vertical 
notions  (mb  hr~')  and  kinematic  isentropic  trajoctorias,  Storm 
position  as  of  1915  GMT, 
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Over   the   2000-2130  GMT  period   the   308K  and    312K  isentropic   sur- 
faces  vertical  motion   patterns    present   sinking  motion  over   most   of 
the   grid,  with    larger  values    (25-30  mb  hr      )    in   the    northern  half   of 
the  grid.      The   320K  and   332K  theta   surfaces   vertical  motions   still 
show  rising    in  front  and    into   the   storm  system,   but   the    patterns   are 
not  as  well  defined.      The   2000-2130  GMT  convergence    fields    on  the    320K 

isentropic  surface   are  well   correlated  with   the   storm  system  and   show 

-4  ~1 

strong  convergence   of   1.0  to    1.5  x   10        sec        in  the   system    (Figs.    4a, 

4b). 

Darkow  (1968)  has  examined  the  total  energy  environment  of  severe 

storms  on  the  synoptic  scale  of  motion  and  has  determined  a  "Total 

Energy  Index".   This  index  is  computed  as  the  difference  between  the 

"total  energy"  (c  T  +  gz  +  Lq)  at  500  mb  and  the  total  energy  at  850 
P 

mb.   Index  values  more  negative  than  -2.0  cal  gm  indicate  severe  storms 
and  associated  tornado  activity  are  highly  probable. 

The  NSSL  isentropic  plane  section  computer  output  for  the  nine 
radiosonde  stations  gives  the  static  energy.   This  data  was  investi- 
gated in  terms  of  a  "Total  Energy  Index".   Values  at  1700  GMT 
indicated  a  region  of  -2.0  in  the  area  where  the  thunderstorms 
began  to  occur.   Values  in  the  same  region  were  -4.0  or  larger 
by  2000  GMT  (Fig.  5).   Values  over  the  remaining  grid  area  were 
slightly  negative  or  positive. 

Total  energy,  defined  by  equation  (8),  was  investigated  for 
storm  evolvement.   The  312K  and  320K  isentropic  surfaces  for  the  2000 
and  2130  GMT  analyses  show  a  region  of  maximum  total  energy  in  the  storm 
vicinity,  with  decreasing  values  toward  the  NW  and  SE.   This  region  is 
seen  in  Figure  6.   As  shown  by  Figure  7,  the  contour  pattern  of  total 
energy  closely  resembles  the  specific  humidity  distribution,  which 
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indicates  the  Lq  term  dominates  the  pattern.   Thus  the  relative  humidity- 
analysis  is  a  critical  factor  in  determining  the  total  energy.   The 
relative  humidity  is  not  the  most  reliable  information  from  a  radio- 
sonde sounding,  particularly  at  higher  altitudes.   Since  the  Lq  term 
seems  to  dominate  the  total  energy  contours,  unless  the  relative 
humidity  field  is  adequately  specified,  the  total  energy  computations 
will  not  provide  good  information  with  which  to  study  the  storm.   There 
are  also  inaccuracies  in  the  Montgomery  stream  function  field  which  may 
amplify  the  error  in  total  energy  computations. 

Fankhauser  (1968)  has  demonstrated  the  wind  field  relationship  for 
this  same  May  28,  1967  storm  system.   Using  chaff  trajectories  and 
aircraft  measurements,  the  streamline  and  isotach  pattern  of  the  com- 
posite relative  wind  field  with  respect  to  the  storm  echo  are  presented 
(Fig.  8).   The  320K  isentropic  surface  almost  corresponds  with  the  500 
mb  isobaric  surface,  thus  a  good  comparison  can  be  made.   In  Figure  9 
the  isentropic  trajectories  show  that  flow  is  being  deviated  around  the 
storm,  convergence  into  the  storm  and  divergence  behind  the  storm.   The 
acceleration  field  (Fig.  10)  also  shows  maxima  to  the  right  and  left  of 
the  storm  and  minima  in  front  of  and  behind  the  storm.   Fankhauser  has 
demonstrated  the  flow  for  an  isolated  storm  echo  and  the  fields  are 
instantaneous.   The  isentropic  trajectories  show  average  conditions 
for  the  larger  storm  system,  and  yet  the  results  are  strikingly  similar. 
The  comparison  times  are  also  different,  which  may  indicate  this  wind 
distribution  is  present  throughout  the  storm  cycle  for  both  isolated 
cells  and  the  entire  storm  complex.   Thus  the  1.5  hour  kinematic  isen- 
tropic trajectories  correlate  very  well  with  the  chaff  trajectories 
that  Fankhauser  has  obtained. 
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Figure  8,   2145  GMT  aimplified  atreainline  and  iaotach  pattern 
relating  the  compoaite  relative  wind  field  with  reapect  to 
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D.   INERTIAL-GRAVITY  WAVE  ASSOCIATION  TO  THE  SEVERE  STORM 

A  close  look  at  various  analyses  seems  to  imply  an  oscillation 
present  in  the  network  area.   Thus  an  investigation  into  the  possibility 
of  a  wave  like  phenomenon  present  over  the  data  area  was  considered  as 
a  possible  source  of  the  storm  activity  and  storm  movement. 

An  oscillatory  tendency  can  be  observed  on  the  320K  theta  surface 
composite  trajectories  (Fig.  11).   Measurements  in  the  center  of  the 
grid  region  show  a  wave  length  of  about  150  KM  and  period  of  about  3 
hours.   Radar  presentations  of  the  most  intense  storm  activity  also 
indicate  a  wave  length  of  about  150  KM  (Fig.  12). 

Surface  pressure  tendencies  over  the  network  area  were  studied  to 
determine  if  the  oscillation  could  be  detected  at  the  surface.  Micro- 
barograph  traces  were  available  at  35  stations  over  the  data  area.   The 
pressure  tendencies  indicated  falling  pressures  of  about  0.03  inches  of 
Hg  per  1.5  hours  over  the  grid  during  the  entire  life  cycle  of  the  storm. 
Figure  13  illustrates  this  tendency  over  the  area.  A  noticeable  fluctu- 
ation was  not  apparent  at  most  stations .   Even  the  presence  of  the  storm 
was  difficult  to  detect  (see  stations  3A,  3C,  k~B   and  kC   on  Fig.  13). 

Synoptic  pressure  tendency  patterns  were  obtained  for  1800  and  2^00 
GMT,  and  a  mean  synoptic  pressure  tendency  over  the  grid  area  was  cal- 
culated.  This  value  was  then  subtracted  from  the  network  stations 
pressure  tendency  at  each  time  frame  in  order  to  calculate  a  "meso- 
scale  pressure  tendency".   The  results,  as  seen  in  Figure  lU,  indicate 
both  rise  and  fall  areas.   The  rise  and  fall  areas  show  a  half  wave 
length  of  about  100  KM  along  the  mean  low  level  wind  field,  and  a  150 
KM  wave  length  between  two  negative  centers  along  a  line  of  mean  wind 
direction. 
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Figure  12,   2045  GMT  storm  poaitioa, 
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Time  sections  at  selected  stations  along  the  mean  wind  direction 
(see  Fig.  15)  were  made  in  order  to  study  the  time  variations  of 
vertical  motion,  static  stability,  relative  humidity,  wind  speed  and 
wind  direction.   The  time  sections  at  stations  LTS,  FSI,  CHK  and  TIK 
(Figs.  l6a-l6d)  indicate  an  apparent  periodicity  in  the  vertical  motion 
and  static  stability  patterns.   The  vertical  motions  are  the  average 
omega  obtained  by  the  trajectory  method,  and  the  static  stability  is 
the  value  centered  between  two  adjacent  theta  surfaces.   Measurements 
on  these  time  sections ,  from  the  vertical  motion  patterns ,  give  a  wave 

length  of  about  150  KM,  wave  period  of  3  hours  and  wave  speed  of  lU  m 

-1 
sec 

It  can  be  seen  in  Figure  l6a  that  a  region  of  vertical  motion  of 
+30  mb  hr~  is  occurring  at  station  LTS  at  1915  GMT,  between  the  308K 
and  312K  theta  surfaces.  With  the  measured  wave  speed  of  ik   m  sec 
(50  KM  hr   ),  a  propagating  wave  without  change  in  amplitude,  constant 
wave  length  and  moving  in  the  direction  of  the  mean  wind  would  produce 
a  region  of  about  +30  mb  hr""""  at  station  FSI  about  20l*5  GMT  and  station 
CHK  about  2lU5  GMT.   This  condition  is  verified  at  these  stations  for 
the  specified  times  (Figs.  l6b  &  l6c). 

With  the  assumption  of  a  wave-like  phenomenon  present ,  a  study  was 
conducted  to  determine  if  such  a  wave  phenomenon  and  other  supporting 
dynamic  and  thermodynamic  forces  could  be  correlated  with  the  onset  of 
the  storm.   The  earliest  radar  presentation  at  1817  GMT  indicated  small 
cumulus  activity  in  the  western  portion  of  the  grid  area.   Time  sections 
at  stations  LTS,  FSI,  CHK  and  TIK  (Figs.  17a-17b)  show  that  at  1700  GMT 
there  is  a  pronounced  shift  in  wind  direction  between  the  stations  on 
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Figure  17a,   Tim*  section  at  station  LTS,   Thin  solid  limes  ere 
thets  surface  positions,   Heary  solid  limes  are  relatire  humidity 
(%),   Dashed  lines  represent  wind  speed  (m  sec"*),   Short  solid 
lines  represent  wind  direction  (360  degs  at  top  of  page), 


& 


1700 


1830  10OO 

Time   thaun) 


2130 


2500 
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the  332K  theta  surface.   This  wind  shift  could  be  further  evidence  of 
the  passage  of  a  wave  through  the  area  around  1700  GMT. 

The  relationship  between  the  horizontal  curvature  of  the  trajectory 
(ICJ  and  the  streamline  (K  )  can  be  related  by  Blaton's  equation: 

hm*.*r*n  (13) 

s 

3  f  f* 
where  V  is  in  the  wind  speed  and  —r—  is  the  local  rate  of  turning  of  the 
s  ot 

3  f  f 
horizontal  wind  direction.   At  station  FSI  (Fig.  1Tb)  — — -  is  positive, 

ot 

thus  K_^  should  be  greater  than  K  ,  which  is  what  is  observed  (Fig.  18). 

3  -pf 
At  station  CHK  (Fig.  17c)  — —  is  negative,  thus  K^  should  be  less  than 

K  ,  which  is  also  the  case  (Fig.  18).   If  the  streamlines  are  wave- 
shaped  and  if  the  wave  moves  with  a  speed  greater  than  the  mean  wind 
speed,  the  curvature  of  the  streamlines  and  trajectories  may  be  of 
opposite  sign  and  the  trajectories  will  have  a  curvature  greater  than 
the  stream! ines.   This  approach  implies  that  a  wave  having  a  phase 
speed  greater  than  the  mean  wind  speed  is  indeed  present  in  the  region. 

The  1700-1830  GMT  vertical  motion,  1700  GMT  static  stability  and 
divergence  patterns  were  superimposed  upon  the  1700  GMT  streamline  flow 
on  the  332K  theta  surface  (Fig.  18).   A  parcel  of  air,  with  respect  to 
a  wave  that  is  traveling  faster  than  the  mean  wind  speed,  is  actually 
moving  toward  the  west.   Thus  the  air  parcels  are  experiencing  con- 
vergence, sinking  motions  and  entering  a  stability  minimum  area  in  the 
western  portion  of  the  grid. 

From  the  conservation  of  potential  vorticity  on  an  isentropic 

surface,  which  is  given  by: 

d_ 

dt  UQa)  =  0  (lU) 
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where  £  and  o  are  defined  as  in  equations  (3)  and  (7),  an  air  parcel 
entering  a  trough  region  would  normally  experience  an  increase  in 
vorticity,  thus  requiring  a  decrease  in  static  stability.   From  the 
continuity  equation  in  isentropic  coordinates  in  the  form: 

f-°Vv  =  0  (15) 

it  is  seen  that,  where  a  is  positive,  a  parcel  experiencing  convergence 
must  have  a  decrease  in  static  stability. 

There  thus  appears  good  evidence  of  the  presence  of  a  wave 
phenomenon  in  the  region  in  the  early  stages  of  storm  activity.   There 
also  is  a  good  indication  that  the  wave  has  a  phase  speed  greater  than 
the  mean  wind  speed.   The  1700  GMT  model  on  the  332K  isentropic  surface 
gives  evidence  that  the  parameters  can  only  be  consistent  for  a  wave 
speed  greater  than  the  mean  wind  speed.   With  this  dynamically  consistent 
model  there  has  been  established  a  region  conducive  to  convective 
activity  for  the  1700-1830  GMT  period  which  could  initiate  the  severe 
storm  activity. 

For  adiabatic ,  frictionless  flow  along  a  streamline,  the  equation 
of  motion  in  isentropic  coordinates  reduces  to  the  form: 

-5F--H  (l6> 

By  assuming: 

where  C  is  the  phase  speed,  equations  16  and  17  can  be  reduced  to  the 

form: 

9M 

C  =  V  +  ||  (18) 
s    9V 
s_ 

3S 
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where  V  is  the  mean  wind  speed.   Now  by  assuming  a  phase  speed  of 
twice  the  mean  wind  speed,  a  relationship  for  adjusting  the  mass 
field  to  the  wind  field  can  be  obtained  in  the  form: 

I  -  vs  sr  (19) 

such  that  AM  3  AV  along  the  streamline  for  the  same  As.   Figure  19 
illustrates  the  adjustment  that  is  necessary  in  the  Montgomery  stream 
function  field.   This  adjusted  field  is  dynamically  consistent  with 
the  model  that  has  been  established  on  the  332K  theta  surface  at  1700 
GMT.   Air  parcels,  moving  with  respect  to  the  wave,  which  are  experi- 
encing an  acceleration  are  also  experiencing  a  decrease  in  Montgomery 
stream  function.   Thus  knowledge  of  the  wave  in  the  region  can  be 
helpful  in  adjusting  the  mass  and  wind  fields  such  that  they  are 
dynamically  consistent. 
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5.   PROBLEMS  ENCOUNTERED 

One  of  the  most  demanding  and  important  parts  of  the  research  are 
the  initial  analyzed  fields.  The  problem  has  required  the  analysis  of 
100  data  fields.   Proper  analysis  involves  both  spatial  and  time  con- 
siderations ,  as  well  as  dynamical  and  thermodynamic al  relationships 
between  variables.  The  horizontal  gradients  are  very  important  as 
interpolations  are  made  between  network  stations.   Subjective  smoothing 
can  destroy  the  very  detail  that  may  well  be  important ,  but  at  the 
same  time  under  smoothing  may  produce  anomalies  that  invalidate  sub- 
sequent studies.   Station  bias  must  be  carefully  considered. 

On  the  synoptic  scale ,  the  normal  concern  of  changes  in  a  certain 
parameter  that  must  be  detected  over  a  specified  distance  or  time  will 
usually  be  expressed  in,  say,  whole  numbers.   However,  for  mesoscale 
work  one  must  detect  changes  that  amount  to  tenths  or  hundredths  for  a 
specified  distance  or  time.   The  instrumentation  for  data  collecting  is 
the  same  for  both  scales  of  motion,  thus  the  sensitivity  of  the  instru- 
ments and  the  magnitude  of  instrument  error  and  measurements  is  very 
critical. 

Engelmann  and  Davis  (1968)  have  found  that  the  standard  observational 

•7 
error  in  radiosonde  data  will  give  a  computational  error  of  0.0U  x  10 

2    -2 
cm  sec   in  the  Montgomery  stream  function.  These  errors  result  from 

errors  in  the  instruments  and  receiving  equipment,  baseline  check  and 

in  the  mean  temperature.   The  M  fields  were  analyzed  for  each  0.05  x 

7   2-2 

10  cm  sec  "  interval,  which  is  about  the  value  of  the  computational 

error.   Thus  the  M  values  and  analyses  fields  contain  inaccuracies 
which  effect  calculations  using  these  fields.  It  was  seen  in  section 
UD  that  the  Montgomery  stream  function  fields  required  rather  large 
adjustments  in  order  to  be  dynamically  consistent  with  the  winds. 
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Figures  20a-20e  indicate  the  changes  in  the  analyses  for  the 
320K  isentropic  surface  at  time  2000  GMT  when  balloon  displacement, 
location  of  the  storm  at  analysis  time,  and  a  conscientious  effort 
toward  spatial  and  time  continuity  were  made.   The  320K  surface  was 
selected  as  a  comparison  for  the  re-analysis  since  its  height  fluctu- 
ations were  minimal  over  the  entire  grid  and  between  each  data  time. 
Isolines  of  Montgomery  stream  function  were  oriented  along  streamlines 
and  rapid  changes  perpendicular  to  the  flow  were  minimized.   Pressure 
isolines  were  also  oriented  along  the  flow,  as  were  relative  humidity 
values.   Also  larger  values  of  relative  humidity  were  analyzed  for, 
in  the  vicinity  of  the  storm.   Wind  direction  and  speed  considerations 
were  carefully  matched,  particularly  with  respect  to  time  continuity. 

The  effects  of  the  re-analyses  on  the  trajectories  and  values  of 
the  associated  parameters  were  investigated  to  determine  if  any  sig- 
nificant changes  resulted.   Figure  21  shows  the  comparison  between 
the  initial  analysis  (IA)  and  the  re-analysis  (FA)  for  selected  tra- 
jectories and  their  associated  parameters.  Vertical  motion  (VM) 
(mb  hr   ),  acceleration  (A)  (m  sec   ),  vorticity  change  (V  CHG)  (sec   ), 
divergence  change  (D  CHG)  (sec   ),  total  energy  change  (TE  CHG) 
(10  erg  gm   ),  static  stability  change  (SS  CHG)  (  K  mb~  )  and  potential 
vorticity  change  (PV  CHG)  (10     K  sec  cm  gm   )  are  given  for  each  of 
the  trajectories.  Vertical  motions  and  accelerations  are  smaller  for 
the  re-analyzed  case  and  there  are  two  cases  which  have  sign  changes. 
The  re-analyzed  trajectories  have  more  cyclonic  curvature  in  the  vicinity 
of  the  storm,  and  trajectory  1  has  more  anticyclonic  curvature.   Other 
parameters  show  little  difference. 
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Using  equation  16  and  the  2000-2130  GMT  accelerations  that  were 
computed  from  the  kinematic  isentropic  trajectories,  the  expected 
changes  in  the  Montgomery  stream  function  values  along  a  selected 
AS  of  50  KM  were  computed  at  various  areas  over  the  grid  on  the  320K 
theta  surface  (Fig.  22).   The  solid  black  lines  are  the  analyzed  M 

field  based  on  the  data  at  the  nine  NSSL  upper  air  stations  at  2000 

7   2-2 
GMT.   The  values  are  given  in  units  of  10  cm  sec   .   The  dashed 

black  lines  represent  the  2000-2130  GMT  accelerations  (m  sec"  hr~  ) 

and  the  short  green  lines  are  the  position  of  AS  over  the  area.   In 

Table  1  a  comparison  of  the  predicted  value  and  measured  value  for 

the  M  values  are  given.   In  many  cases  there  is  almost  an  order  of 

magnitude  difference,  with  the  predicted  values  having  a  much  lower 

absolute  value.  Using  these  predicted  changes  and  approximate  values 

at  each  station,  an  adjustment  to  the  M  field  was  made  (Fig.  22  -  red 

dashed  lines).   This  new  M  field  is  now  dynamically  consistent  with 

the  winds.   It  would  thus  seem  that  the  Montgomery  stream  function 

and  hence  the  temperature  and  height  fields  are  in  error.   It  is  not 

conclusive  however,  but  along  with  the  adjustments  in  the  M  field 

discussed  in  Section  kD   there  appears  to  be  unacceptable  errors  in 

the  initial  data  fields  of  Montgomery  stream  function.   From  this 

point  of  view,  therefore  it  appears  the  wind  data  is  more  accurate. 
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TABLE  1 
Predicted  and  Measured  M  Value  Comparison 


M 

M 

Calculation 

Acceleration 

Predicted 

Measured 

Point 

(m  sec  ""  hr   ) 
-l.U 

(      2    ~2\ 

(cm  sec   ; 

+  .020 

(      2    -2\ 
(cm  sec   ) 

1 

-.08 
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-0.9 

+  .013 

-.07 

3 
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+  .001 
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-0.U 
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+1.1 
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+  .01U 
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+  .05 

IT 

-0.9 

+  .013 

+  .01 
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-1.5 

+  .021 

-.ok 

19 

-0.7 

+  .010 

-.03 

20 

+1.0 

-.01^ 

0.0 
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6.   COMPARISON  WITH  THEORY 

In  the  basic  assumption  of  frictionless ,  adiabatic  flow  along 
an  isentropic  surface,  on  which  the  study  was  based,  there  are 
three  parameters  which  should  be  conserved  for  an  air  parcel  moving 
in  the  flow.   These  parameters  are  potential  temperature  (0)  (  K), 
specific  humidity  (q)  (gm  gm   )  and  potential  vorticity  ( f* )  (10~    K 
sec  cm  gm   ) .   The  initial  input  fields  into  the  computer  program  were 
analyzed  on  constant  theta  surfaces,  thus  the  program  constrained  theta 
to  be  conserved  for  each  trajectory.   However,  specific  humidity  and 
potential  vorticity  were  allowed  to  vary. 

Selected  composite  trajectories  were  studied  in  order  to  determine 
both  short  (1.5  hour)  and  long  (6  hour)  term  variations.   In  Figure  22 
and  Table  2  the  deviations  are  shown  for  the  320K  theta  surface. 
Vertical  motions  (u?)  (mb  hr   ),  accelerations  (A)  (m  sec  "  hr   ), 
static  stability  (<T)  (  K  mb   )  and  total  energy  (TE)  (10  erg  gm"  ) 
were  also  studied  in  conjunction  with  the  conserved  properties  to 
determine  if  any  correlations  exist.   Trajectories  labeled  1A-1F  are 
for  1700-1830  GMT;  2A-2F  for  1830-2000  GMT;  3A-3F  for  2000-2130  GMT; 
and  1+A-l+F  for  2130-2300  GMT.   Both  initial  (i)  and  final  (f)  values 
for  each  of  the  properties  are  given. 

With  the  exception  of  1A-1F  trajectories,  the  specific  humidity 
was  somewhat  conserved  for  each  trajectory.   Changes  are  random  and 
do  not  appear  to  be  associated  with  vertical  motions  or  accelerations. 
Trajectories  2E  and  3F  indicate  their  deviations  may  be  from  the  storm 
effects,  since  their  initial  point  is  outside  the  storm  and  final  point 
inside  the  storm  region. 
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TABLE  2   (Continued) 
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7.   CONCLUSIONS 

The  research  has  demonstrated  the  critical  nature  of  mesoscale  data 
both  in  instrumentation  and  analysis  errors.   It  was  found  that  the  wind 
data  were  perhaps  the  most  reliable  for  use  in  the  study.   Subsequently 
the  1.5  hour  kinematic  isentropic  trajectory  relationship  to  the  severe 
local  storm  warranted  considerable  exploration.   It  was  further  noted 
that  the  320K  theta  surface  (approximately  500  mb)  showed  the  tra- 
jectory and  air  parcel  properties  which  were  more  closely  associated 
with  the  life  cycle  of  the  storm. 

The  trajectories  indicated  that  prior  to  storm  development  the  wind 
veered  with  height  at  all  levels ,  flow  was  anticyclonic  with  the  ex- 
ception of  cyclonic  flow  in  the  NW  corner  in  upper  levels,  and  a  wind 
maximum  region  existed  in  the  NW  corner  of  the  network.   With  the 
presence  of  an  inertial-gravity  wave  at  the  332K  theta  surface  at  1700 
GMT,  the  dynamically  consistent  vertical  motions,  static  stability, 
convergence,  trajectories,  streamlines  and  adjusted  Montgomery  stream 
function  fields  support  conditions  which  could  initiate  convective 
activity.   It  was  also  noted  that  the  storm  is  to  the  NE  of  the 
Wichita  mountains,  and  the  existance  of  a  lee  wave  supporting  the 
inertial-gravity  wave  is  worthy  of  additional  studies. 

The  study  of  various  air  parcel  properties  on  the  isentropic 
surfaces  has  provided  insights  into  the  storm  dynamics.   The  1830- 
2000  GMT  vertical  motions  indicated  supporting  rising  motions  for 
storm  growth.   The  2000-2130  GMT,  320K  theta  surface  convergence 
patterns  help  to  illustrate  the  intensity  of  the  storm  during  this 
period.   The  lack  of  supporting  vertical  motions,  and  the  anti- 
cyclonic  flow  in  each  level  tended  to  provide  the  response  necessary 
for  decay  of  the  storm  over  the  2130-2300  GMT  period. 
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The  orientation  of  the  storm  along  the  upper  level  trajectories 
in  the  early  stages  and  along  the  low  level  trajectories  in  the  latter 
stages  may  be  indicative  of  storm  movement.   When  the  storm  orientation 
is  considered  along  with  the  inertial-gravity  wave,  there  exists  the 
possibility  the  storm  formed  at  high  levels  and  propagated  downward 
into  the  lower  levels.   It  is  concluded  that  the  kinematic  isentropic 
trajectories,  over  short  periods  of  time,  do  provide  an  effective 
method  for  studying  local  severe  storm  activity  in  the  mesoscale  range. 
8.   SUGGESTIONS  FOR  FUTURE  DEVELOPMENTS  AND  IMPROVEMENTS 

Since  the  initial  data  fields  are  critical  to  the  success  of  the 
computer  model,  it  would  be  of  benefit  to  investigate  the  application 
of  error  analysis  techniques  in  order  to  improve  the  initial  data. 
Further  studies  are  necessary  when  corrections  for  Montgomery  stream 
function  changes  can  be  made,  in  order  to  improve  the  final  trajectories. 
Incorporation  of  diabatic  effects  in  regions  where  air  parcels  are  be- 
coming saturated  prior  to  the  final  trajectory  point  can  be  included  in 
the  model.   Further  studies  of  different  storm  systems  is  encouraged  in 
order  to  test  the  short  period  kinematic  isentropic  trajectories,  and 
to  determine  if  they  can  be  utilized  as  a  forecasting  aid.   Since  there 
existed  deviations  in  conservation  of  specific  humidity  and  potential 
vorticity  along  the  trajectories,  additional  research  will  be  necessary 
to  investigate  the  conservation  of  these  properties.   With  the  sup- 
porting evidence  of  an  inertial-gravity  wave  in  the  region,  a  theoreti- 
cal approach  is  recommended  in  order  to  study  the  phenomenon  in  the 
range  of  wave  length  and  wave  speed  encountered  in  this  case  study. 
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APPEUDIX  A  -  ISOTROPIC  KINEMATIC  TRAJECTORIES 

This  appendix  includes  the  kinematic  trajectories  which  have  been 
computed  for  the  308K,  312K,  320K  and  332K  isentropic  surfaces  for 
the  time  periods  1700-1830,  1830-2000,  2000-2130  and  2130-2300  GMT. 
They  were  produced  by  the  objective  computer  model  for  the  May  28, 
1967  severe  storm  case.   The  grid  is  an  x,  y  grid,  with  equal  grid 
spacing. 
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APPENDIX  B  -  ISENTROPIC  TRAJECTORY  COMPUTER  PROGRAM 

This  appendix  includes  the  objective  computer  model.   The  program 
is  written  in  FORTRAN  language.   Each  routine  contains  comment  cards 
which  explain  the  routine  and  each  step. 
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